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Summary 

Electrostatic  field  effects  in  high  pressure  dis¬ 
charges  which  provide  physical  insight  into  the  pro¬ 
blem  of  electrode  design  are  derived.  The  coupled  set 
of  charge  transport  equations  and  Poisson's  equation 
have  been  analytically  studied  using  perturbation 
theory  techniques.  Our  study  reveals  that  the  homoge- 
eous  field  must  be  supplemented  by  the  inhomogeneous 
contributions  of  space  charge  and  polarization  effects 
which  arise  from  nonuniform  plasma  ionization  and 
local  nonconservation  of  each  species  current.  The 
validity  of  a  commonly  used  boundary  condition  is 
examined  and  found  untenable.  Implications  for  the 
computational  design  of  electrodes  are  presented. 

Introduction 

The  rational  engineering  design  of  electrodes  for 
externally  ionized  as  well  as  self-sustained  glow 
discharges  requires  an  appreciation  of  dominant 
physical  effects  as  well  as  mathematical  consequences 
of  the  modelling  equations.  Contemporary  design 
practice  invokes  techniques  which  vary  in 
sophistication  from  the  widely  used  empirical  method 
of  cut  and  try  to  numerical  methods  of  solving  the 
nonlinear  Poisson  equation  in  arbitrary  materials  and 
geometries . 

To  some  investigators,  a  very  appealing  feature  of 
discharge  physics  is  its  insatiable  appetite  for 
computationally  recondite  theories  which  must  of 
course  be  self-consistently  linked  together.  Unfor¬ 
tunately,  in  our  haste  to  compute,  we  sometimes  lose 
sight  of  the  most  elementary  issues  which  can  then 
sabotage  us  from  undergrowth  of  arcanna  in  which  they 
hide.  Such  a  problem  in  electrode  design  and  in 
calculations  of  the  discharge  field  is  how  to  pre¬ 
scribe  the  boundary  conditions  of  the  field.  In 
tackling  this  problem,  which  is  addressed  below,  we 
have  begun  with  a  fresh  look  at  the  electrostatics  of 
discharges  based  on  elementary  undergraduate  consider¬ 
ations.  Our  conclusions  have  provided  us  with  intui¬ 
tive  insight  into  the  physical  mechanisms  which  govern 
the  field  behavior  and  have  guided  us  in  the  construc¬ 
tion  of  our  computational  models  used  in  electrode 
design. 

In  the  following,  we  limit  ourselves  to  the  design  of 
two  dimensional  electrode  profiles.  In  practice,  an 
approximate  solution  to  the  real  problem  of  deter¬ 
mining  a  three  dimensional  profile  may  be  obtained  by 
smoothly  fitting  together  rectilinear  and  axi- 
symmetric  profiles.  Theoretical  work  is  in  progress 
to  address  the  computational  problem  of  a  self- 
consistent  three  dimensional  electrode  design. 


Historical  Survey 

The  Zeoroeth  Order  and  most  prevalent  method  of 
contemporary  electrode  design  assumes  that  the  field 
is  completely  homogeneous  and  a  resort  is  made  to  the 
conformal  mapping  solutions  of  Laplace's  equation 
which  appear  in  Maxwell's  treatise.  The  algorithm  for 
design  was  worked  out  by  Rogowski  and  appears  in 
Cobine's  book.  The  Rogowski  profile  results  from  a 
semi- inf inite  cartesian  coordinate  system.  The 
correct  solution  for  a  finite  width  coutour  is  in 
principle  available  but  contains  implicit  functions  of 
the  elliptic  functions  of  the  first  and  second  kind. 
The  Chang  profile  results  from  analytically  continuing 
the  Rogowski  profile  and  requiring  maximal  uniformity 
of  the  field  on  the  center  of  the  elctrode.  In  addi¬ 
tion  to  these  analytic  techniques,  rather  elaborate 
computer  codes  now  exist  to  solve  Laplace's  equations 
in  arbitrary  materials  and  subject  to  standard  mixed 
boundary  conditions  on  a  finite  enclosing  surface. 

The  objective  of  these  analytical  techniques  is  to 
produce  an  electrode  profile  on  which  the  field 
intensity  is  maximally  uniform  and  monotonically 
decreasing  from  the  symmetry  plane  of  the  profile. 
However  theoretical  and  experimental  evidence  strongly 
suggests  that  this  may  not  be  the  most  desirable 
engineering  criterion.  In  ref  k  investigations  are 
reported  in  which  the  electrode  profile  was  shaped  to 
provide  for  abrupt  deviations  of  the  field  from  the 
above  conditions  near  the  edge  of  the  discharge.  The 
result  was  a  reduction  of  discharge  constriction,  a 
very  sharp  discharge  edge  and  very  flat  optical 
uniformity.  Furthermore  it  is  a  commonly  observed 
fact  that  large  discharge  devices  tend  to  arc  in  pre¬ 
ferential  locations  on  the  electrodes.  These  loca¬ 
tions  are  usually  found  where  one  can  intuitively 
perceive  the  existence  of  large  field  stress  concen¬ 
trations.  Thus  in  such  cases  device  failure  is  due  to 
a  local  weakness  in  design  and  not  a  global 
deficiency.  A  highly  desirable  engineering  design 
objective  is  then  to  distribute  such  local  stresses 
throughout  a  larger  area  to  improve  uniformity  of  per¬ 
formance  throughout  the  device. 

There  are  two  issues  which  are  raised  by  the  above 
techniques.  A  major  limitation  of  the  analytic 
techniques  and  of  a  standard  "Laplace"  computer  code 
is  that  they  are  normally  limited  In  validity  to 
rectilinear  coordinate  systems  with  translational 
symmetry  in  the  third  dimension  and  thus  may  not  be 
casually  applied  to  cylindrical  geometries  which  are 
essential  to  switch  design.  The  solutions  derived 
from  conformal  mapping  techniques  may  not  be  applied 
to  cylindrical  geometries  because  the  applicable  part 
of  the  cylindrical  Laplacian  is  not  harmonic  in  the 
sense  of  complex  variable  theory.  The  "Laplace"  code 
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solutions  on  the  other  hand  do  not  in  general  use  the 
cylindrical  Laplacian  for  the  numerical  development. 
Thus  an  uncritical  use  of  standard  profiles  in 
cylindrical  devices  which  involve  a  mere  rotation  of  a 
rectilinear  solution  about  the  symmetry  axis  is 
unwarranted.  Work  is  now  in  progress  to  develop 
standard  techniques  for  this  problem. 

The  second  major  limitation  which  comes  to  light  in  a 
practical  application  of  the  numerical  technique  is 
that  a  real  device  geometry  may  occasionally  possess  a 
boundary  which  is  sufficiently  distant  from  the  cavity 
that  it  sets  up  an  asymptotic  boundary  condition.  In 
practice  the  investigator  may  then  impose  an  artifi¬ 
cial  boundary  on  the  device  which  approximates  the 
asymptotic  consequences  in  a  bounded  region.  However, 
even  though  this  artificial  condition  may  in  some 
sense  be  a  good  overall  approximation  to  the  true 
field  behavior,  it  can  introduce  distortions  to  the 
fringing  field  behavior  and  hence  to  the  roll-off 
contour  of  an  electrode  profile.  A  better  practical 
approximation  to  the  problem  of  distant  boundaries  is 
to  impose  the  true  asymptotic  conditions  in  a  finite 
domain  which  has  been  mapped  from  the  "infinite" 
domain.  The  technique  of  numerical  mapping  as  it  is 
addressed  in  the  following  companion  paper  thus 
retains  the  essential  fringing  field  behavior 
characteristics  of  the  conformal  mapping  solutions. 

The  technique  of  merely  solving  Laplace's  equation  for 
electrode  design  works  well  in  self-sustained  dis¬ 
charges.  However,  as  we  will  show  below,  inhomo¬ 
geneous  source  terms  do  not  completely  disappear  for 
this  case,  they  are  only  less  important  than  those  in 
the  externally  sustained  case. 

In  externally  sustained  discharges,  the  spatially  non- 
uniform  primary  ionization  causes  the  cavity  field  to 
deviate  substantially  from  the  homogeneous  field. 
Instead  of  solving  Laplace's  equation,  a  nonlinear 
Poisson  equation,  resulting  from  current  conservation 
requirements,  is  solved  numerically.  In  our 
analysis  we  will  investigate  assumptions  which  enter 
into  this  procedure.  We  will  find  that  the  assumption 
of  vanishing  normal  field  at  the  edge  of  a  discharge 
confinement  region  is  not  merely  invalid  as  would  be 
generally  suspected,  but  that  it  is  also  entirely 
wrong  and  misleading.  Next  we  will  investigate  the 
assumption  of  electron  drift  current  conservation 
which  leads  to  the  applicable  nonlinear  Poisson 
equation.  We  will  discover  nonconservation  of 
electron  drift  current  with  a  consequent  space  charge 
and  second  order  field  distortion. 


Here  n.  (x)  is  the  number  density  of  ions  or 
i ,  e 

electrons,  v.  (E/N)  is  the  drift  velocity  due  to  the 
1 » e 

neutral  gas,  N,  and  the  normalized  electric  field, 

+  ,-4  .  +  .  .  . 

E  (x)  =  -  V  $.  The  recombination  rate  is  a(E/N)  and 
the  combined  secondary  ionization  and  attachment  rate 
is  3(E/N).  The  spatially  varying  primary  ionization 
source  term,  S(x)  ,is  also  an  implicit  function 
of  E(x)  and  B(x)  .  But  we  neglect  the  self-consistent 
coupling  of  S  to  E  and  B  for  this  analysis.  In  the 
event  that  S  is  provided  for  by  an  electron  beam,  then 
the  fast  primary  electrons  will  eventually  be 
thermalized  by  the  discharge  gas  and  will  contribute 
to  the  space  charge  term.  This  can  be  a  substantial 
effect,  but  we  will  neglect  it  here  for  the  sake  of 
clarity  in  our  following  arguments.  For  the  sake  of 
algebraic  simplicity  we  simplify  our  considerations  to 
the  case  of  3=0.  Analagous  conclusions  may  be  drawn 
from  an  extension  of  our  analysis  to  the  case  of 
3*0.  We  have  neglected  the  diffusion  current  and 
its  derivatives  in  comparison  to  the  convection 
current  and  its  derivatives  for  the  specialized  case 
of  high  pressures. 

The  conventional  procedure  is  to  assume: 

^ni  e 

a)  Steady  state,  — =  0 

b)  Quasi-neutrality,  n_.=  n  ,  for  all  x;  and 

c)  Separate  species  current  conservation, 

4;  4-  ‘  -4  4 

V  •  n  v  =  0  =  V  •  n,  v. 
e  e  i  i 

From  the  electron  current  equation,  for  example,  we 
then  obtain  n^=  /S/a,  which  together  with  current 
conversation  yields  a  nonlinear  conductivity  equation: 


Analysis 


(E,x 


We  will 

derive 

our  higher 

order  effects 

from 

the 

coupled 

set  of 

equations 

for  electron 

and 

ion 

transport  together  with  Poisson's  equation. 


3n 

e  £  4- 

ttt - H  V«n  v  =  S  -  an.n  +  Bn 

3 1  e  l  e  e 


3n . 

14- 

-r— - h  V«n.  v.  =  S  -  an.n  +  Bn 

3 1  li  l  e  e 


2  e 

V  $  = - (n n  ) 

e0  i  e 


This  further  results  in  the  nonlinear  Poisson 
equation : 

V2*  - 

Here  v.  (E/N)  and  a(E/N)  are  known  either 
l ,  e 

experimentally  or  by  kinetic  calculations.  The 
nonlinear  Poisson  equation  is  then  solved  subject  to 
appropriate  boundary  conditions  to  give  a  First-Order 
solution  to  the  discharge  field  problem. 

In  the  case  of  a  self-sustained  discharge,  the 
external  ionizing  term,  S(x) ,  disappears  and  must  be 
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replaced  by  the  combined  secondary  ionization  and 

attachment  term  in  .  We  then  obtain  n.  =  3/a  which 

i,e 

upon  inserting  into  the  current  conservation  equation 
and  applying  the  chain  rule  yields: 


3 

3  E 


3  ( E/N)  v. 

i  ,6 


(E/N) 


E  a  (E/N) 


■}  V-E 


Which  results  in  Laplace's  equation  for  the 
field,  V2$  =  0  . 


E(x) 


f  s(x)1  C 
f  s(x)]  ^ di 


and  thus 


P 


_e  v  ™  r _ 

o  o  l 

I  [s(x)J  ?d£ 


Thus  the  First-Order  development  of  the  coupled  set  of 
equations  for  a  self-sustained  discharge  results  in  a 
homogeneous  electric  field.  Therefore  electrodes  for 
this  case  may  be  based  on  a  harmonic  analysis  as  has 
been  the  practice. 


The  charge  stored  within  a  stream  tube  between  elec¬ 
trodes  is  thus 

sq)]~g-[s(o)]~r-  | 

/2[s(x)]  r’dZ  ) 


In  the  following  we  challenge  an  assumption  which  is 

2-5 

often  and  recklessly  made.  Namely  that  the 

electric  field  normal  to  a  discharge  confinement  wall 
vanishes.  This  assumption  is  equivalent  to  assuming 
that  the  fringing  fields  which  occur  with  any  bounded 
conductor  in  free  space  are  contained  by  the  confined 
plasma  "since  the  plasma  is  a  conductor."  Another 
argument  one  may  occasionally  encounter  which  supports 
the  idea  of  confined  field  lines  in  the  plasma  is  that 
the  plasma  is  so  strongly  polarized  that  it  traps  the 
field  lines  by  dielectric  focusing  as  in  the  optical- 
phenomenon  of  total  internal  reflection.  We  will 
analyze  both  points  of  view  and  demonstrate  that  the 
assertion  of  a  vanishing  normal  component  of  the  field 
is  insupportable  and  spurious. 

Space  Charge  Effects  From  Current  Conservation 

For  convenience,  we  model  a  and  vg  as  power  functions 
of  E/N  and  assume  that  N  is  constant.2  That  is  we 
neglect  gasdynamic  effects. 

Let  a  =  A2E  2a  and  v  =  BE^ 
e 

and  define  s(x)  -  and  5 


Using  the  definition  of  capacitance  per  unit  area  C/A 
=  (Q/A )/V  we  find  that  the  ratio  of  the  capacitance  of 
the  discharge  to  that  of  the  free  space  capacitor  is 
given  by  the  tern  in  braces.  For  a  uniform  source  we 
see  that  the  ratio  vanishes.  For  a  very  severe  case, 
we  may  approximate  s(x)  =  exp(-kx)  .  This  corresponds 
to  a  real  AFWL  device,  PULSAR,  where  k  =  -0.216,  a  = 
0.78,  b  =  0.56.  We  obtain  the  ratio  of  the  capaci¬ 
tance  of  a  flux  tube  to  that  of  the  free  space 
capacitor  plates  to  be  /C  =  p  =  0.081.  That  is 
the  free  space  capacitor  contributes  over  90%  of  the 
field  behavior  to  the  device  and  hence  the  field  is 
mostly  homogeneous  and  mostly  fringes  at  its  edges. 


We  pause  to  note,  however,  that  even  though  the  bulk 
effect  of  the  discharge  on  the  overall  capacitance  is 
small,  the  local  field  intensity  may  deviate  substan¬ 
tially  from  the  homogeneous  field.  In  the  case  of 
PULSAR  we  have  (A  =  5) : 


E(x) 

E 

o 


pA  tix  .20 
M«.  ,  e  I  0.80 

e  -1 


Polarization  Effects  for  Current  Conservation 


Then  the  space  charge  density  is: 

p  =  £  ?*E  =  — £  ?  E-V  ln[  s(x)  ]  . 

o  o 

Now  from  i  =  i  =  en  v  =  e/S  rr  Ea+^ 
o  e  e  e  A 

l 

and  V  =  -  J  E  •  d£  ,  we  obtain 
o 

o 


Since  we  now  know  that  the  discharge  stores  charge  in 
the  bulk  we  would  like  to  know  how  it  behaves  in  small 
localities.  In  a  differential  volume  of  cross  sec¬ 
tional  area  dA  and  length  dx  along  a  flux  line  we  have 
a  charge  dQ  and  a  voltage  drop  dV  =  Edx.  Together 
these  lead  to  an  incremental  capacitance  per  unit  area 
of 
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c' 


dQ/dA 

dV 


p  dx 


E  dx 


E*7  In  [s(x)] 


accomplished  by  the  mapping  techniques  presented  in 
the  following  companion  paper. 


We  thus  see  that  our  externally  ionized  discharge  can 
store  charge  and  possesses  an  incremental  capacitance 
due  to  its  non-uniform  conductivity.  We  will  now 
outline  a  calculation  of  the  discharge's  associated 
equivalent  relative  permittivity  on  the  basis  of 
electrostatic  energy  arguments.  We  have 


A  second  consequence  of  the  smallness  of  the  distor¬ 
tions  introduced  by  the  polarization  of  the  discharge 
is  that  we  must  now  include  all  surrounding 
dielectrics  in  our  calculations. 

Nonconservation  of  Species  Currents 


U  -  l/2e  E2  =  l/2e  E2 

O  0  0 


t[s(x)l  ?  ~|2 
i 

-/  [s(x>]  ?d2- 


1/2  e  _e  E 
rel  o  o 


The  term  in  braces  is  thus  an  equivalent  relative 
local  dielectric  constant  of  the  discharge  in  the 
sense  of  indicating  either  constriction  or  dilation  of 
the  field  lines  in  the  transverse  direction  as 
Faraday's  rubber  lines  of  force  would  tend  towards. 
For  the  case  of  PULSAR, 
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rel 


(x) 


r  Ml  iixi  2  'f  1 .50 

'  y2  ,  e  M  0.65 
e  -1 


n- 


Although  we  are  presenting  an  exaggerated  and  pathol¬ 
ogical  case  we  can  clearly  see  the  general  trends. 
The  discharge  acts  like  a  dielectric  lens  which  bends 
the  fringing  field  lines  away  from  the  cavity  near  the 
cathode  and  towards  the  cavity  near  the  anode.  This 
explains  the  empirically  determined  engineering 
practice  of  sharply  rounding  cathode  roll-off  contours 
and  mildly  rounding  anode  roll-off  contours  when  the 
e-beam  is  injected  from  the  cathode  side. 


In  this  section  we  outline  an  investigation  into  the 
validity  of  assuming  conservation  of  each  species 
current.  Unfortunately  the  lengthy  and  tedious 
algebra  involved  mandates  a  condensation  of  our 
reasoning  into  its  most  salient  features.  We  will 
find  that  each  species  currant  is  not  independently 
conserved  and  that  as  a  result  the  true  field  will 
vary  from  the  first  order  field  by  Second-Order 
effects  due  to  power  loading. 

First  we  note  that  in  deriving  our  non-linear  Poisson 
equation,  that  we  could  just  as  easily  have  used 
vi  as  ve*  However  the  dependence  of  ve  on  E/n  varies 
widely  among  gases  with  a  power  dependence  varying 
from  b  =  0.5  for  CO2  to  b  =  0.9  for  N2  whereas  for 
most  bases  the  power  for  is  very  close  to  unity. 
We  are  thus  led  to  seek  out  disparities  in  the  overall 
behavior  due  to  this  difference. 


b  c 

Let  v  =  BE  and  v .=  CE 
e  i 


then 


£  Fc-b 

B  E 


Once  again  we  see  that  the  discharge  only  slightly 

affects  the  true  field  in  comparison  to  the  homo¬ 
geneous  field  (although  significantly  for  electrode 

design  purposes)  .  Thus  in  the  vicinity  of  the  dis¬ 
charge  edge  we  are  completely  unjustified  in  asserting 
that  the  normal  electric  field  component  vanishes  as 
is  done  in  ref  5-7.  The  effect  of  imposing  this 
artificial  condition  is  to  distort  the  curvature  of 
the  potential  surfaces  near  the  roll-off  thus  inval¬ 
idating  any  electrode  profile  design  accomplished  in 
this  manner.  In  order  to  retain  the  correct  field 

behavior  near  the  roll-off  we  must  correctly  match  the 
field  behavior  in  the  asymptotic  region.  This  is  best 


where 


For  the  pathological  case  of  PULSAR  we  have  S  <  1CT-3 

and  S'  =  k  J  <  10~3  .  We  may  thus  neglect  6  in  our 
subsequent  remarks,  although  In  general  &  must  be 
retained  throughout  all  manipulations  to  ensure  that 
correct  perturbation  orders  are  addressed. 

From  jj,  +  je  =  jQ,  we  obtain 


280 


j 


o 


C  Fc“b 

B  E 


■  - 1 


1  +  f 


In  the  spirit  of  successive  approximations  we  thus  see 

that  V  •  j.  =0  implies  j  V*E  |  >  0  which  in  the 
l  ,e 

next  higher  order  treatment  implies  that 
I  V.  j.  J  >  0. 


Using  E  in  terms  of  j^/jQ  from  above,  we  may  rewrite 

a(E)  and  v.  and  thus  n.  in  terms  of  l./i  • 

v  '  i,e  i,e  Ji  Jo 

Thus  the  ion  current  conservation  equation  may  be  cast 
into  a  form  which  depends  only  on  j.  /j  and 

auxilliary  parameters: 


where 


n=2 


a+b 

c-b 


and  0 


2 


a+c 

c-b* 


We  have  analytically  investigated  this  equation  with 
perturbation  theoretic  techniques  based  on  the 
assumption  of  a  slowly  varying  j^(x)  and  using 
linearization  wherever  possible  due  to  the  smallness 

of 


permittivity .  Thus  this  power  loading  induced  space 
charge  effect  cannot  in  general  be  ignored. 


Having  established  the  effect  of  species  current  non¬ 
conservation  on  the  discharge  we  have  now  justified 
the  need  to  solve  the  fully  coupled  set  of  transport 
equations  together  with  Poisson's  equation  in  order  to 
correctly  determine  the  electric  field  for  electrode 
design  for  externally  ionized  discharges.  This  aspect 
of  the  general  problem  is  currently  being  pursued  with 
a  mid-term  priority. 


Finally,  we  may  note  that  even  if  we  set  S(x)=0  in  our 
rewritten  ion  conservation  equation,  that  j^/j  will 
still  have  a  spatial  dependence  (we  of  course  must  add 
the  g  term  which  we  dropped  earlier) .  Thus  to  second 
order,  even  self-sustained  discharges  will  possess 
electric  fields  which  deviate  from  the  homogeneous 
solution. 


Conclusion 


The  electric  field  of  an  externally  ionized  discharge 
is  mostly  homogeneous.  Space  charge  can  affect  the 
local  field  intensity  by  as  much  as  20  percent  in 
severe  cases.  Discharge  polarization  can  enhance  or 
diminish  transverse  field  structure,  i.e.,  fringing, 
by  as  much  as  50  percent.  A  rational  electrode  design 
theory  must  predict  these  effects. 
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